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Abstract 

This papeiQ deals with low maximum likelihood (ML) decoding complexity, full-rate and full-diversity 
space-time block codes (STBCs), which also offer large coding gain, for the 2 transmit antenna, 2 receive 
antenna (2 x 2) and the 4 transmit antenna, 2 receive antenna (4 x 2) MIMO systems. Presently, the best 
known STBC for the 2 x 2 system is the Golden code and that for the 4 x 2 system is the DjABBA 
code. Following the approach by Biglieri, Hong and Viterbo, a new STBC is presented in this paper for 
the 2x2 system. This code matches the Golden code in performance and ML-decoding complexity for 
square QAM constellations while it has lower ML-decoding complexity with the same performance for 
non-rectangular QAM constellations. This code is also shown to be information-lossless and diversity- 
multiplexing gain (DMG) tradeoff optimal. This design procedure is then extended to the 4x2 system 
and a code, which outperforms the DjABBA code for QAM constellations with lower ML-decoding 
complexity, is presented. So far, the Golden code has been reported to have an ML-decoding complexity 
of the order of M"^ for square QAM of size M. In this paper, a scheme that reduces its ML-decoding 
complexity to M^\/M is presented. 

L Introduction And Background 

Multiple-input multiple-output (MIMO) transmission has attracted a lot of interest in the last decade, 
chiefly because of the enhanced capacity it provides compared with that provided by the single-input, 
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single-output (SISO) system. The Alamouti code [1] for two transmit antennas, due to its orthogonality 
property, allows a low complexity ML-decoder. This scheme led to the development of the generalized 
orthogonal designs [2]. These designs are famous for the simplified ML-decoding that they provide. 
They allow all the symbols to be decoupled from one another and hence, are said to be single-symbol 
decodable. Another bright aspect about these codes is that they have full transmit diversity for arbitrary 
complex constellation. However, the limiting factor of these designs is the low code rate (refer Section 
HI] for a definition of code rate) that they support. 

At the other extreme are the well known codes from division algebra, first introduced in [3]. The well 
known perfect codes [4] have also been evolved from division algebra with large coding gains. These 
codes have full transmit diversity and have the advantage of a very high symbol rate, equal to that of the 
VBLAST scheme, which, incidentally doesn't have full transmit diversity. But unfortunately, the codes 
from division algebra including perfect codes have a very high ML-decoding complexity (refer Section 
HI] for a definition of ML-decoding complexity), making their use prohibitive in practice. 

The class of single-symbol decodable codes also includes the codes constructed using co-ordinate 
interleaving, called co-ordinate interleaved orthogonal designs (CIODs) [5], and the Clifford-Unitary 
Weight single-symbol decodable designs (CUW-SSD) [6]. These designs allow a symbol rate higher than 
that of the orthogonal designs, although not as much as that provided by the codes from division algebra. 
The disadvantage with these codes when compared with the Orthogonal designs is that they have full 
transmit diversity for only specific complex constellations. 

The Golden code [7], developed from division algebra, is a full-rate (see Section JI] for the definition 
of full-rate), full-diversity 2x2 code for integer lattice constellations, but has been known to have a high 
ML-decoding complexity, of the order of M^, where M is the size of the constellation used (it is shown in 
Section rvll] that this can be reduced significantly to M^a/M when the constellation employed is a square 
QAM). It has to be mentioned that the codes presented in [8] and [9] also have the same coding gain and 
ML-decoding complexity as the Golden code does. With a view of reducing the ML-decoding complexity, 
two new full-rate, full-diversity codes for QAM constellations have been proposed for the 2x2 MIMO 
system. The first code was independently discovered by Hottinen, Tirkkonen and Wichman [10] and by 
Paredes, Gershman and Alkhansari [11], which we call the HTW-PGA code and the second, which we 
call the Sezginer-Sari code, was reported in [12] by Sezginer and Sari. Both these codes enable simplified 
ML-decoding (see Section JI] for a definition of simplified ML-decoding), achieving a complexity of the 
order of in general, and for square QAM (shown in Section I VIII ). These codes have a slightly 
lower coding gain than the Golden code and hence show a slight loss in performance compared to the 
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Golden code. These codes sacrifice the coding gain for simplified ML-decoding complexity. 

For 4 transmit antennas, the popular codes are the quasi-orthogonal designs, first introduced in [13] 
and the CIOD for 4 transmit antenna [5], both of which are rate one codes. The CIOD is known to 
be single symbol decodable and the MDC-QOD [14] is also single symbol decodable. But when 2 or 
more receive antennas are employed, these codes cannot be considered to be full-rate. The perfect code 
for 4 transmit antennas has a high rate of 4 complex symbols per channel use but its use in practice 
is hampered by its high decoding complexity, even with the use of sphere decoding [15], [16]. For a 
4x2 MIMO system, the best performing code has been the DjABBA code [10], which beats even the 
punctured perfect code for 4 transmit antennas in performance [17], [18]. This code was designed for 
performance alone and has a high ML-decoding complexity, of the order of Af^, as shown in Section 
IVIII The first attempt at reducing the ML-decoding complexity for a 4 x 2 system while maintaining 
full-rate was made by Biglieri, Hong and Viterbo [17]. The full-rate code that they have proposed, which 
we call the BHV code, has an ML-decoding complexity of the order of for general constellations, 
(though this has been reported to be A/^ in [17]), but doesn't have full-diversity. However, the code 
matches the DjABBA code in the low SNR scenario and betters the punctured perfect code in codeword 
error performance (CER). 

The contributions of this paper are as follows 

• we propose a new full-rate, full-diversity STBC for the 2x2 MIMO system. This code has an 
ML-decoding complexity of the order of Af'^ in general, as compared to A/^ for the Golden code. 
For square QAM, the ML-decoding complexity of our code is of the order of A'/^\/2, the same as 
that of the Golden code. 

• Our code also matches the Golden code in coding gain for QAM constellations and is shown to 
have the non-vanishing determinant (NVD) property for QAM constellations and hence, is DMG 
optimal. We also show that our code is information-lossless. 

• We propose a new full-rate, full-diversity STBC for 4 x 2 MIMO systems, having ML-decoding 
complexity of the order of for arbitrary complex constellations, and of the order of M^^/M for 
square QAM constellations, whereas the corresponding complexity for the DjABBA code are 
and respectively. It also has a higher coding gain than the DjABBA code for 4- and 16-QAM 
constellations and hence, a better CER performance. 

• We state the conditions that allow simplified ML-decoding and show that for square QAM constel- 
lations, the ML-decoding complexity of the Golden code can be reduced to Af^\/M. 
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The remaining content of the paper is organized as follows : In Section |lll we give the system model 
and the code design criteria. In Section JIIJ we present our code for the 2x2 MIMO system and show 
that it is information-lossless. In Section |IVl we show that our code has the NVD property and DMG 
optimality. In Section |Vl we present our code for the 4x2 MIMO system. Section |VI] deals with the 
low complexity ML-decoding of these codes. In Section IVIII we analyze the ML-decoding complexity 
for the Golden code, the HTW-PGA code, the DjABBA code and the BHV code. The simulations results 
constitute Section IVIIII Concluding remarks are made in Section |lXl 

Notations: Throughout, bold, lowercase letters are used to denote vectors and bold, uppercase letters 
are used to denote matrices. Let X be a complex matrix. Then X^, X^ and det [X] denote the transpose, 
Hermitian and determinant of X, respectively. For a complex variable x, xj and xq denote the real 
and imaginary part of x, respectively. Also, j represents and the sets of all integers, all real and 

complex numbers are denoted by Z, M and C, respectively. The Frobenius norm and the trace operations 
are denoted by \\.\\f and tr [.] respectively. The operation of stacking the columns of X one below the 
other is denoted by vec(K). The Kronecker product is denoted by 0, It and Ot denote the TxT identity 
matrix and the null matrix, respectively. The inner product of two vectors x and y is denoted by (x, y). 
For a complex random variable X, X J\fc{0,N) denotes that X has a complex normal distribution 
with mean and variance A^. For any real number m, rnd[m] denotes the operation that rounds off m 
to the nearest integer, i.e.. 



rnd[m] 



[m\ if [m] — m > m — [m\ 
[m] otherwise 



For a complex variable x, the (.) operator acting on x is defined as follows 



- A 
X = 



XI -XQ 
XQ XI 



The (.) can similarly be applied to any matrix X G ([^nxm j-gpiacing each entry Xij by Xij, i 



1, 2, • • • , n, j = 1, 2, • • • , m , resulting in a matrix denoted by X G 
Given a complex vector x = [xi,X2, - ■ ■ , x is defined as 



p 2n X 2 m 



and x' is defined as 



X = [XII,XIQ, ■■■ , Xnl, Xuq]'^- 



x' = [-XiQ,Xu, , • • • , -XnQ, Xnlf^- 
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It follows that X = [x x']. 

II. System Model 

We consider Rayleigh quasi-static flat-fading MIMO channel with full channel state information (CSI) 
at the receiver but not at the transmitter. For nt x n,. MIMO transmission, we have 

Y = HS + N (1) 

where S € C"*^^ is the codeword matrix, transmitted over T channel uses, N G (C">-xT ^ complex 
white Gaussian noise matrix with i.i.d entries ~ A/c (0, A^o) and H G C""^"' is the channel matrix with 
the entries assumed to be i.i.d circularly symmetric Gaussian random variables Ac (0, 1). Y G ([^"■r-xT 
is the received matrix. 

Definition 1: (Code rate) If there are k independent complex information symbols in the codeword 
which are transmitted over T channel uses, then, the code rate is defined to be k/T complex symbols 
per channel use. For instance, for the Alamouti code, k = 2 and T = 2. So, its code rate is 1 complex 
symbol per channel use. 

Definition 2: (Full-rate code). An STBC is said to be fiill-rate if it transmits at the rate of rimin 
complex symbols per channel use, where rirain = min {nt,nr). 

So, the Alamouti code can be considered to be full-rate for 2 x 1 MIMO alone, while the Golden code 
is full-rate for n^, > 2. 

Considering ML-decoding, the decoding metric that is to be minimized over all possible values of 
codewords S is given by 

M(S) = ||Y-HS||| (2) 

Definition 3: (Decoding complexity) The ML decoding complexity is a measure of the maximum 
number of symbols that need to be jointly decoded in minimizing the ML decoding metric. This number 
can be k in the worst scenario, k being the total number of information symbols in the code. Such a 
code is said to have a high ML-decoding complexity, of the order of M^, where M is the size of the 
signal constellation. If the code has an ML-decoding complexity of order less than M^, the code is said 
to admit simplified ML-decoding. For some codes, all the symbols can be independently decoded. Such 
codes are said to be single-symbol decodable. 

Definition 4: (Generator matrix) For any STBC S that encodes k information symbols, the generator 
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matrix G is defined by the following equation [17] 

vec{S) = Gs. (3) 

where s = [si, 52, • ■ ■ ) -Sfc]^ is the information symbol vector. 

An STBC can be expressed in terms of its weight matrices (linear dispersion matrices) as follows 

k 

S = ^ A2i-.iSii + A2iSiQ (4) 

i=l 

Here, Aj, i = 1, 2, • • • ,2k are the weight matrices for S. It follows that 

G = [uec(Ai) t)ec(A2) ••• i;ec(A2fc)] (5) 

It is well known [21], that an analysis of the PEP leads to the following design criteria: 

1) Rank criterion: To achieve maximum diversity, the codeword difference matrix (S — S) must have 
full-rank for all possible codeword pairs and the diversity gain is ntUr- If full-rank is not achievable, 
then, the diversity gain is given by rrir, where r is the minimum rank of the codeword difference 
matrix over all possible codeword pairs. 

2) Determinant criterion : For a full ranked STBC, the minimum determinant 5min, defined as 

H' 



A_ 

■'mm 



s -s s -s 



(6) 



should be maximized. The coding gain is given by (5mm)^^"% with nt being the number of transmit 
antennas. 

If the STBC is non full-diversity and r is the minimum rank of the codeword difference matrix over 
all possible codeword pairs, then, the coding gain 6 is given by 



S = min Aj 



s-s 



\i=l 



where Aj, i = 1, 2, • • • , r, are the non-zero eigen values of the matrix — — . It should be 

noted that for high signal-to-noise ratio (SNR) values at each receive antenna, the dominant parameter 
is the diversity gain which defines the slope of the CER curve. This implies that it is important to first 
ensure full-diversity of the STBC and then try to maximize the coding gain. 
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III. The Proposed STBC For 2x2 MIMO And Information-Losslessness 

In this section, we present our STBC [22], [23] for 2 x 2 MIMO system. The design is based on the 
CIODs, which were studied in [5] in connection with a general class of single-symbol decodable codes 
which includes complex orthogonal designs as a proper subclass. Specifically, for 2 transmit antennas, 
the CIOD is as follows. 

Definition 5: The CIOD for 2 transmit antennas [5] is 



X(si,S2) 



(V) 



SlI + jS2Q 

52/ + jSlQ 

where G C, i = 1, 2 are the information symbols and sn and s-iq are the in-phase (real) and quadrature- 
phase (imaginary) components of Si, respectively. Notice that in order to make the above STBC full- 
rank, the signal constellation A from which the symbols Sj are chosen should be such that the real part 
(imaginary part, resp.) of any signal point in A is not equal to the real part (imaginary part, resp.) of any 
other signal point in A [5]. So if QAM constellations are chosen, they have to be rotated. The optimum 
angle of rotation has been found in [5] to be ^tan~^2 radians and this maximizes the diversity and 
coding gain. We denote this angle by 9g. 
The proposed 2x2 STBC S is given by 



S(xi,X2,X3,X4) =X(si,S2) + e-'^X (53 , S4) P 



(8) 



where 



The four symbols si, 82,83 and 84 G A, where A is a. 6 g radians rotated version of an integer QAM 
signal set, denoted by Aq, which is a finite subset of the integer lattice, and xi,X2, X3, X4 G Aq, i.e. 



Xi ; 



1,2,3,4. 



1 

1 



• P is a permutation matrix designed to make the STBC full-rate and is given by P = 

• The choice of 9 in the above expression should be such that the diversity and coding gain are 
maximized. We choose 9 to be 7r/4 and show in the next section that this angle maximizes the 
coding gain. 

Explicitly, our code matrix is 



S(xi,X2,X3,X4) 



SlI + jS2Q eJ''/^(s3/ + jS4Q) 



e^'"/^(s47+jS3Q) 



S2I + jSlQ 



(9) 
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with Sii = cos{0g)xu - sin{6g)xiQ and siq = sin{9g)xii + cos{9g)xiQ, i = 1,2,3,4. 

The minimum determinant for our code when the symbols are chosen from the regular QAM constel- 
lations (one in which the difference between any two signal points is a multiple of 2) is 3.2, the same 
as that for the Golden code, which is proved in the next section. The generator matrix for our STBC (as 
defined in Definition 01), corresponding to the information vector consisting of symbols follows: 



G 



C0S{9g) 


—sin{9g) 


























sin{9g) 


COS{9g) 


























siniOg) 
V2 


COS (6*3 ) 


COS (6*3 ) 


sin{8g) 
V2 














sin{9g) 
V2 


COS{0g) 


COS{0g) 

V2 


sin{8g) 
V2 














COS{9g) 


sin(9g) 
V2 


sin(9g) 
V2 


COSffig) 














cos(eg) 


sin(9g) 
V2 


sin{8g) 
V2 


COS{9g) 








COS{9g) 


—sin{9g) 














sin{9g) 


cos{9g) 





















(10) 



It is easy to see that this generator matrix is orthonormal. In [9], it was shown that a sufficient condition 
for an STBC to be information-lossless is that its generator matrix should be unitary. Hence, our STBC 
has the information-losslessness property. 



IV. NVD Property And DMG Optimality Of The 2x2 Code 

In this section we show that the proposed code has the NVD property [7], which, in conjunction with 
fuU-rateness, means that our code is DMG tradeoff optimal [19]. We also show that the angle 7r/4 in ^ 
maximizes the coding gain. 

Theorem 1: The minimum determinant of the proposed 2x2 code, given by when the symbols 
are chosen from is 1/5. 

Proof: The determinant of the codeword matrix S can be written as 



det{S) = {sii + js2Q){s2i + jsio) - j[{szi + jsAQ){sAi + jszq)]. 
Using Sii = {si + s*)/2 and jsjg = (sj - s*)/2 in equation we get. 



(11) 



4def(S) = (si + Si + S2 - S2)(S2 + S2 + Sl — S*i) — j[(s3 + S3 + Si - S4)(S4 + S4 + S3 — S3)] 

= ((Sl + S2) + (Si - S2)*) ((Si + S2) - (Si - S2)*) - i[((s3 + S4) + (S3 - S4)*) ((S3 + S4) - (S3 - S4)*)]. 
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Since Si = e^^^Xi,i = 1,2,3,4, with Si ^ A, Xi Aq, a subset of defining A = {xi + X2), 
B ^{xi- X2)*, C = (xs + X4) and D = (x^ - X4)*, with A, B, C and D e we get 

Metis) = [e^^-A + e^^'^^B){e^^^ A - eT^^-B) - j[{e^^^C + e-^'^^ D){e^^^C - e^^^'W)] 

Since e^'^^s = cos{2eg) + jsin{29g) = (1 + 2j)/^/b, we get 

iVbdetiS) = (1 + 2j)iA^ - jC^) - (1 - 2j){B^ - jD^). (12) 

For the determinant of S to be 0, we must have 

{l + 2j){A^-jC^) = {l-2j){B^-jD^) 
^{l + 23)\A^-3C^) = 5iB'-jD^). 

The above can be written as 

Al-jCf = 5{B^-jD^) (13) 

where Ai = {1 + 2j)A, Ci = (1 + 2j)C and clearly Ai,Ci £ Z[j]. It has been shown in [8] that (fT3] ) 
holds only when Ai = B = Ci = D = 0, i.e., only when xi = X2 = x^ = X4 = 0. This means that 
the determinant of the codeword difference matrix is only when the codeword difference matrix is 
itself the zero matrix. So, for any distinct pair of codewords, the codeword difference matrix is always 
full-rank for any constellation which is a subset of Also, the minimum value of the modulus of 
the R.H.S of (ini) can be seen to be 4. This occurs for {A,B,C,D) = (1,1,0,0) or (0,0,1,1). The 
occurrence of any other combination of A,B,C and D that results in a lower value of the modulus of 
the R.H.S of (fT2l) can be ruled out after noting that xi,X2,X3 and X4 take only values from For 
eg. {A, B,C, D) = (l,j, 0, 0) is one such combination, but it is easy to see mathematically that such a 
combination cannot occur for Xi e = 1,2,3,4. So, \det{S)\ > l/Vb, meaning that the minimum 

determinant for the code is 1/5. ■ 
In particular, when the constellation chosen is the regular QAM constellation, the difference between 
any two signal points is a multiple of 2. Hence, for such constellations, |(iet(S-S')| > 4/\/5, where S 
and S' are distinct codewords. The minimum determinant is consequently 16/5 and hence the proposed 
code has the NVD property [7]. Now, from [19], where it was shown that full-rate codes which satisfy 
the NVD property achieve the optimal DMG tradeoff, our proposed STBC is DMG tradeoff optimal. 
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As a byproduct of Theorem [T] we arrive at the following lemma. 

Lemma 1: The choice of 7r/4 for 6* in (HJl maximizes the coding gain of the proposed 2x2 code for 
QAM constellations. 

Proof: Consider the CIOD whose codeword has the structure shown in (|7]). The set of codeword 
difference matrices of the CIOD is a subset of the set of the codeword difference matrices of the proposed 
2x2 code, whose codeword structure is given in Q. It is to be noted that the minimum determinant 
and hence the coding gain of a code depend on the codeword difference matrices of the code. In (fT2l) . if 
we let C = L> = 0, we arrive at the expression for the determinant of a codeword matrix of the CIOD. 
So, for the CIOD, whose codeword matrix is denoted by S, we have 

4.Vhdet{S) = (1 + 2j)A^ - (1 - 2j)s2_ (24) 

where, A = (xi +2:2) and B = {xi — X2)* , with xi and X2 taking values from It is evident that the 
minimum of the modulus of the R.H.S of ([141 ) is 4, which occurs for ^ = = 1. So, the minimum of 
the absolute value of the determinant of a codeword matrix of the CIOD when the symbols take values 
from (not all taking zero values) is 1/a/5. When the symbols take values from the regular QAM 
constellation, the minimum of the absolute value of determinant of a non-zero codeword difference matrix 
is 4/\/5 and hence, the minimum determinant for the CIOD is 16/5. We have already shown that the 
minimum determinant for our 2x code is 16/5, when the symbols take values from the regular QAM. 
This shows that the choice of 7r/4 for in dSjl indeed maximizes the coding gain. ■ 

V. The Proposed STBC For The 4x2 MIMO System 

In this section, we present our STBC for the 4x2 MIMO system [24] following the same approach 
that we took to design the 2x2 code. The design is based on the CIOD for 4 antennas, whose structure 
is as defined below. 

Definition 6: CIOD for 4 transmit antennas [5] is as follows: 

SlI+jS'iQ -S2I+jSAQ 

S2I + jS4Q Sil - jS3Q 

S37 + jSiQ -SaI + JS2Q 

S47 + jS2Q S23 - jSiQ 

where Si,i = 1, • • • ,4 are the information symbols as defined in the previous section. Here again, the 
symbols are chosen from a rotated version of the regular QAM constellation, with 9g being the angle of 



X(si, 52,53,54) 



(15) 
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rotation. 

The proposed STBC is obtained as follows. Our 4x4 code matrix, denoted by S, encodes eight symbols 
xi, • • • ,xs drawn from a QAM constellation, denoted by Aq. As before, we denote the rotated version 
of Aq by A. Let Sj = e^^<'Xi, z = 1, 2, • • • 8, so that the symbols Si are drawn from the constellation A. 



The codeword matrix is defined as 



S{xi,X2, ■■■ ,X8) = X(si,S2, 53,54) + e-'^X(s5,S6,S7,S8)P 



(16) 



with 6 e [0, 7r/2] and P being a permutation matrix designed to make the STBC full-rate and given by 

10 
1 
10 
10 

The choice of 9 is to maximize the diversity and coding gain. Here again, we take to be 7r/4. This value 
of 9 provides the largest coding gain achievable for this family of codes. This is so because the minimum 
determinant for the CIOD as defined in ([TS] ) (which can also be obtained by letting the variables 55, 
sq, s-j and ss be zeros in ([T6l )) is 10.24 [14] for unnormalized QAM constellations. The value of the 
minimum determinant for our 4x2 code, obtained for unnormalized 4-QAM and 16-QAM constellations 
is 10.24, which was checked by exhaustive search. This shows that the choice of 7r/4 maximizes the 
coding gain. The resulting code matrix is as shown below. 



sii + js3Q -S21 + js^Q e^'^/^is^i + Jsyq) 

S21 + js^Q su - js3Q e^''/^(s67 + Jssq) 

e^^'^^issi + Jsqq) e^''/^(s7/ - jssg) 54/ + js2Q 



-Sil + jS2Q 
S3I - jSlQ 



This code is full-rate only for the 4x2 MIMO system, unlike the perfect space time code [4], which 
is full-rate for 7ir > 4. Also, the generator matrix for our code can be checked to be non-unitary. So, our 
STBC for 4 X 2 MIMO system is not information-lossless. 



VI. Low Complexity ML-Decoding Of The 2x2 And 4x2 Codes 

In this section, we show how our codes admit simplified ML-decoding. The information symbols are 
assumed to take values from QAM constellations. In the general setting, it can be shown that ^ can be 
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written as 

t^^^) = UegX + wec(N) (17) 

where Heg G ]R2n,Tx2fc given by 

Heg=(lT0H)G (18) 
with G G i^2ntTx2A: ^^gjj^g )-j,g generator matrix as in Definition HI so that vec{S) = Gx. and 

X= [XII,XIQ,--- ,XkI,XkQ]'^ 

with Xi,i = 1, • • • ,k drawn from Aq, which is the regular QAM constellation. Using this equivalent 
model, the ML decoding metric can be written as 

M (x) = Wvlic^) - Hegxf (19) 

On obtaining the QR decomposition of Hgg, we get Hgg = QR, where Q G ]^2n^Tx2A: orthonormal 
matrix and R G M^'^^^^' is an upper triangular matrix. The ML decoding metric now can be written as 

M(x) = ||Q^v15y) - Rxf = ||y' - Rxf (20) 

where y' — [y'l, - ■ ■ , y2k] — Q^vec{Y). If Heg = [hi h2 • • • h2k], where hi,i = 1,2, ■ ■ ■ ,2k are column 
vectors, then Q and R have the general form obtained by Gram — Schmidt process as shown below 

Q = [qi q2 Qa • • • q2fc] 

where qj, i = 1, 2, • • • ,2k are column vectors, and 

R^ 



where ri = hi, = = hj - T,]=\{^j^^i)^j^ = i = 2,3,--- ,2k. 

Lemma 2: Let M = [f i f 2 • • • fn][gi g2 • • • ?,nV , where fj = [fn, fi2, ■■■ , fin], gi = [9ii,9i2, ■■■ , 9in] 
eW'<\i = l,2,--- ,n. Then, tr(M) = Er=i(g.,fi>- 



rill (qi,h2) (qi,h3) ... (qi,h2fc) 

||r2|| (q2,h3) ... (q2,h2fe) 

llrall ... (q3,h2fc) 

... Ilrsfcll 
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Proof: From the definition of the trace operation, we have 

n n 

trim) = Y.Y.f^l3^, 

n n n 
i=l j=l i=l 

Theorem 2: For an STBC with k independent complex symbols and 2k weight matrices A;,/ = 
1,2, •• • ,2k, if, for any i and j, i ^ j,l < i,j < 2k, AiAj^ + AjA^ = 0„j, then, the i*^ and the 
j*^ columns of the equivalent channel matrix H^q are orthogonal. 

Proof: We note that the following identities hold for matrices A e C""^",B G €"^^^,0 G C^^" 
and vectors x G C^^Sz G C^^^ 

A = BC4»A = BC (21) 
(z,x) = (z',x') (22) 

With these identities, we proceed as follows 

AiAf + AjAf = On, ^ UAiAfU^ + HA^Af H" = 0„,. 

Applying the (.) operator and using (|2TI) . we get 

HAi (HA, f + HA,- {HAif = 02„, (23) 

This indicates that the real matrix M = HAj(HAj)^ is a skew-symmetric matrix and hence its diagonal 
elements are zeros. Let Aj = [aj^i ai^2 ••• ^i,T], where ai^k,k = 1, 2, • • • , T are the columns of Aj. 
Then, Aj = [a^^i a.[i • • • a^^r y]. Therefore, 

M = [Ha,,i Ha',,1 • • • Ha,,r Ha^,T][Ha,-,i Ha;. ^ • • • Ha,,T m'^^^f 
Since M is real and skew-symmetric, tr{M) = 0. So, 

T 

^{(Ha,,^,Ha,,„,) + (Ha;.^,Ha:,J} = (24) 

m=l 

T 

^ 2Y,{^^3,mMi,n.) = (25) 

m=l 
T 

.-. ^(Haj-^,Ha,,™) = (26) 

m=l 
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where, (l24l) follows from Lemma |2] and (l25l ) follows from (l22l ). Now, 



H 



eg 



(It H) G 



(It H) 



Ha 



1,T 



ai,i a2,i 

ai,2 ^2,2 



a2fc,i 
a2fc,2 



ai,r a2,T • • • ^2k,T 

Hai,i Ha2,i ••• Ha2fc,i 
Hai 2 Ha2,2 • • • Ha2fc,2 



Ha 



2,T 



Ha 



2fc,T 



o 



From the above structure, it is readily seen that for any i j,l < i, j < 2k, if AjA^ +AjA 
then the z*^ and the j*^ columns of Hgg are orthogonal. This follows from ( [26l ). ■ 
Now, let us consider the proposed STBC for 2 x 2 MIMO system. Here, k = A,T = 2. It can be 
verified that the following holds true for l,m G {1, 2, 3, 4} 



A^Af + A/ A, 



H 



o 



(27) 



V/ / m, m + 1, if m is odd 
\/l m,m — 1, if m is even. 

To be precise, dl?]) holds for G {(1, 3), (1, 4), (2, 3), (2, 4)}. Therefore, from Theorem |2l (hi, ha) = 

(hi,h4) = (h2,h3) = (h2,h4) =0. 

Using the above results in the definition of the R-matrix, it can easily be shown that (q^jha) = 
(qi,h4) = (q2,h3) = (q2,h4) = 0. So, the structure of the R-matrix for our 2x2 code is as follows. 



R 



rill 


(qi,h2) 








(Qi^hs) 


(Qiihe) 




(qi,h8 





111*211 








(q2.h5) 


(q2,h6) 


(q2ih7) 


(q2.h8 








Iksll 


(q3.h4) 



















l|j*4|| 


(q4.h5) 






(q4.h8 














Iksll 























Ikell 
























Ikrll 

























Iksll 



(28) 
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The structure of the R-matrix enables one to achieve simpUfied ML-decoding. This is because once 
the symbols X3 and X4 are given, xi and X2 can be decoded independently. In the ML-decoding metric, 
it can be observed that the real and imaginary parts of symbol xi are entangled with one another but 
are independent of the real and imaginary parts of X2 when 2:3 and X4 are conditionally given. So, the 
number of metric computations required is at most (M^)(2M) = 2M^ and hence, the ML-decoding 
complexity is of the order of M^. When the constellation employed is a square QAM so that the real 
and the imaginary parts of each symbol can be decoded independently, the ML-decoding complexity 
can be further reduced as follows. Let x = • • • ,X4q] denote the decoded information vector. 

Assuming that sphere decoding is employed (sphere decoding can be employed for constellations like 
square or rectangular QAM and not for any arbitrary constellation which is a finite subset of the 
following strategy is employed - 

1) A 4 dimensional real SD is done to decode the symbols X4 and x^, and there are such pairs 
for an M-QAM constellation. 

2) Next, X2Q is decoded in parallel with xig, and there are \/M possibilities for each of them. 
Following this, xu and X21 are decoded using hard-limiting, as follows 

xi/ = mi?i|max(2rn(i[-^] , — M) , m| (29) 

X21 = min^max (2rnd^-^^ , — , (30) 

where, 

4 

■'^i - {y'l - r{i,2)XiQ - ^{r{i,2i-i)Xii + ?-(i,2i)^*Q))A(i,i) 

1=3 

4 

U2 = {y's - r(3A)S:2Q - X](^(3,2i-l)^i/ + ?'(3,2i)^iQ))/^(3,3) 

1=3 

where, for simplicity, we have denoted the {i,jY^ entry of the R-matrix by r(^ijy 
So, the ML-decoding complexity of our code for square QAM is of the order of M^\/M- If, however, 
the QAM constellation used is not a square QAM, and cannot be represented as the Cartesian product 
of two PAM constellations (like the 32-QAM constellation, the optimum representation of which is as 
shown in Figure |5]l, then the method described above cannot be employed. So, in such a scenario, the 
ML-decoding complexity becomes M'^, because one requires to decode wholly the complex symbols xi 
and X2, when x-^ and X4 are given. 
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Now, let us consider the proposed STBC for 4 x 2 MIMO system. For this case, = 8, T = 4. It can 
be verified that the condition in (l27l) holds true for l,m G {1,2, •• • ,8}. Hence, from Theorem [2l for 
/, m G {1, 2, • • • , 8} , we have 



(hz,h^) = < 



V/ 7^ m, m + 1, if m is odd 
V/ 7^ m, m — 1, if m is even. 

Using the above result, it can be easily be verified that for Z, m G {1, 2, • • • , 8}, / < m, 

V/ / m, m + 1, if m is odd 
V/ / m, m — 1, if m is even. 

For simplicity, let us define the R matrix as follows 



R 



Ri R2 
Og R3 



where, Ri,R2 and R3 G 



^ then, Ri can be seen to have the following structure 



rill 


(qi>h2) 























111*211 


























Iksll 


(q3'h4) 























l|l*4|| 


























Iksll 

























Ikell 


























llrrll 

























Iksl 



(31) 



The above structure of the matrix R allows our code to achieve simplified ML-decoding as follows - having 
fixed the symbols x^,xq, xj and xg, the symbols xi,X2, x-^ and X4 can be decoded independently. In the 
decoding metric, it can be observed that the real and imaginary parts of symbol xi are entangled with one 
another but are independent of the real and imaginary parts of X2, X3 and X4 when xs, xe, X7 and xg are 
conditionally given. Similarly, X2, x^ and X4 are decoupled from one another although their own real and 
imaginary parts are coupled with one another. So, in general, the ML-decoding complexity of our code is 
of the order of M^. That is due to the fact that jointly decoding the symbols x^,xq, x-j and xg followed 
by independently decoding xi,X2,X3 and X4 in parallel requires a total of (M^)(4M) = 4M^ metric 
computations. However, when square QAM is employed, the ML-decoding complexity can be further 
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reduced as follows. Let x = [xij,xiq, ■ ■ ■ ,x^q\ denote the decoded information vector. Assuming the 
use of a sphere decoder, 

1) an 8 dimensional real SD is done to decode the symbols x^,XQ^x^ and xg. 

2) Next, xiQ, X2Q, X3Q and x^q are decoded in parallel. Following this, xu, X21 x^i and x^j are 
decoded using hard limiting as follows 

Xu = min^max{^2rnd[^^ , — M) , m| (32) 

X21 = min^max[2rnd[-^], —M) , (33) 

i|maa;(2rn(i[y],-M),M| (34) 

i{max(2rnd[y],-M),M} (35) 



X31 = mm< 



X4^i = mmi 



where, 



8 



Ul - {y'l - r{l,2)XlQ - ^(?^(l,2i-l)^i/ + '^(l,2i)^iQ)) A(l,l) 

8 

U2 = {y's - r(3,4)^2Q - X]('^{3,2i-l)^i/ + ?^(3,2i)iiQ)) A{3,3) 

i=5 

8 

i=5 
8 

i=5 

where, r^j j) denotes the (i, jY^ entry of the R-matrix. So, in all, we need to make a maximum of 4M^\/M 
metric computations only. Hence, for square QAM constellations, the ML-decoding complexity of our 
code is of the order of M^\/M. 

VII. Comparison of Ml-decoding complexity of our codes with known 2x2 and 4x2 

STBCS 

The ML-decoding complexity of our 2x2 code was shown in the previous section to be of the order of 
M^. This was due solely to the behavior of the weight matrices which resulted in the R-matrix structure as 
in (l28l) for our 2x2 code. For any code, the weight matrices entirely define the ML-decoding complexity. 
For eg., all the weight matrices of the Alamouti code satisfy the condition in Theorem |2l and hence, 
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the equivalent channel matrix Hgg is orthogonal. So, the R-matrix for the Alamouti code is a diagonal 
matrix and this results in an ML-decoding complexity of the order of M for general constellations. In 
the special case of the constellation being a square QAM constellation, the real and imaginary parts of 
each symbol can be further decoded independently using hard-limiting and the decoding complexity of 
the Alamouti code for square M-QAM constellations is constant. For the Golden code, the ML-decoding 
complexity has been considered to be of the order of in the literature [11], [12], [17]. However, the 
ML-decoding complexity of the Golden code can be reduced to the order of M^^/M for square QAM 
constellations. It can be easily verified, by studying the weight matrices and using Theorem |2l that the 
Golden code has the following R-matrix structure :- 

aOaOaaaa 
OaOaaaaa 
OOaOaaaa 
OOOaaaaa 
OOOOaaaa 
OOOOOaaa 
OOOOOOaa 
0000000a 

where 'a' denotes a possible non-zero entry. This structure makes the ML-decoding complexity of the 
Golden code evident. In general, the ML-decoding complexity is of the order of M^. However, when 
square M-QAM is employed, the following decoding strategy can be employed, assuming that a sphere 
decoder is used. 

1) A 4-dimensional real SD is done to decode the symbols and X3. 

2) Next, X21 and X2q are decoded in parallel. Following this, xu and xiq are decoded as follows 



R 



Golden code 



Xll 



where, 



ui 



|max(2rnd[y] , -M), m| 
xiQ = min^max(2rnd[-^] , — M) , m| 

4 

{y'l - r{l,3)X2I - Xl(^(l:2i-l)^i/ + ^'(l,2i)^iQ))/?'(l,l) 
i=3 



(36) 
(37) 
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A / / 

U2 = m 



'{2,4 



)X2Q - '^{r{3,2i-l)XiI + ?'(3,2j)^jQ)) A(2,2) 



i=3 



^HTW-PGA 



where, as usual, we have denoted the (i, j)*'* entry of the R-matrix by r(j and x = [xu, xiq, • • • , x^q] 
denotes the decoded information vector. 

So, the ML-decoding complexity is of the order of M^\/M, the same as that for our 2x2 code. However, 
for non-rectangular QAM constellations, the Golden code does not admit simplified ML-decoding. The 
codes presented in [8], [9] and [20] also have their R-matrix structures identical to that of the Golden 
code and hence offer the same ML-decoding complexity. 

Considering the HTW-PGA code, the R-matrix structure is observed to be as foUows:- 

aOOOaaaa 
OaOOaaaa 
OOaOaaaa 
OOOaaaaa 
OOOOaOOO 
OOOOOaOO 
OOOOOOaO 
0000000a 

where 'a' again denotes a possible non-zero entry. From this structure, the order of the ML-decoding 
complexity can be easily calculated for the different QAM constellation types. For square M-QAM, it 
is of the order of Af^. This follows from the fact that when the symbols x-^ and X4 are fixed, xu, 
xiQ, X21 and X2Q can be decoded independently from one another and each of them can be decoded 
by using hard-limiting, hence requiring a total of only 4Af^ computations. For non-rectangular QAM 
constellations, the ML-decoding complexity is of the order of 2M^. The Sezginer-Sari code also has a 
similar ML-decoding complexity. The above observations are all captured in Tables |II1 In the table, the 
ML-decoding complexity given for each code is the maximum number of metric computations needed. 

The ML-decoding complexity of our 4x2 code was shown to be of the order of for general 
constellations, and M^\/M for square QAM constellations. This simplified complexity was facilitated 
by the structure of the R-matrix, a part of which had the structure as in (|3T]) . The ML decoding complexity 
of the DjABBA code is of the order of M'^ in general, and of the order of for square M-QAM. To 
the best of our knowledge, this hasn't been mentioned in literature. To see this, one has to look at the 
R-matrix structure for the DjABBA code which, as mentioned before, is dictated by the weight matrices 
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R 



l,DjABBA 



for the code. The structure of the Ri -matrix for the DjABBA code, one corresponding to (|3TI ). is as 
follows 

aOOOaaaa 
OaOOaaaa 
OOaOaaaa 
OOOaaaaa 
OOOOaaaa 
OOOOOaaa 
OOOOOOaa 
0000000a 

where 'a' corresponds to a possible non-zero entry. For square Af-QAM, it is evident that xu, xiq, X21 
and X2Q can be decoded independently from one another, by using hard-limiting, when the symbols x^, 
X4, X5, xq, xj and are fixed. This allows an ML-decoding complexity of the order of M^, with 
metric computations for decoding the other 6 symbols. This scheme can be employed only for square 
QAM constellations, so that the real and the imaginary parts can be decoded independently. However, 
for non-rectangular QAM constellations, one must decode xi and X2 independently, when the rest of the 
symbols are given. So, the ML-decoding complexity is of the order of M^. 

The BHV code, which was designed primarily for simplified ML-decoding complexity, has a complexity 
of the order of in general and of the order of M'^\fM specifically for square M-QAM (Incidentally, 
the authors of [17] haven't claimed this!). This follows from the structure of the Ri-matrix as shown in 
(l38l) . with a denoting a possible non-zero entry. 

For square M -QAM, the following strategy can be employed to decode the symbols. 

1) An 8-dimensional real SD is employed to decode the symbols x^, xq, x-j and xg. 

2) Following this, 3:3/, x^q, xh and x^q are decoded in parallel. Next, xu, xiq, X21 and X2q are 
decoded by employing hard-limiting. 

Hence, the ML-decoding complexity of the BHV code is of the order of M'^^/M, because a maximum 
of AM'^^/M metric computations need to be done in minimizing the ML-decoding metric. But for non- 
rectangular QAM constellations, the pairs (xi, X3) and (x2, X4) have to be decoded in parallel after jointly 
decoding the last four symbols, thus accounting for an ML-decoding complexity of the order of M^. 
Table |lll] captures the ML-decoding complexities for the three codes for the different classes of QAM 
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constellations. 



R 



l,BHV 
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VIII. Simulation Results 



(38) 



In aU the simulation scenarios in this section, we consider quasi-static Rayleigh flat fading channels 
and the plots are shown for the codeword error rate (CER) as a function of the SNR at each receive 
antenna. 



A. 2 X 2 MIMO 

Figure [T] shows the CER performances of our 2x2 code, the Golden code and the HTW-PGA code, 
with all the codes employing the 4 QAM constellation. Figure |2] shows the CER plots for the three codes, 
with the constellation used being 16 QAM. In both the plots, we see that the CER curve for our 2x2 
code is indistinguishable from that of the Golden code and this is due to the identical coding gains of 
the two codes. The HTW-PGA code has a slightly worse performance because of its lower coding gain. 
Table U gives a comparison between the minimum determinants of some well known 2x2 codes. It 
is to be noted that in obtaining the minimum determinants for these codes, we have ensured that the 
average energy per codeword is uniform across all codes, but the average energy per constellation has 
been allowed to increase with constellation size, or in other words, the average constellation energies 
haven't been normalized to unity. 

B. 4x2 MIMO 

Figure [3] shows the CER performance plots for our 4x2 code, the well known DjABBA code [10] 
and the BHV code [17], with all the codes using the 4-QAM constellation. Figure |4] shows the CER 
performance for 16 QAM. Both the plots exhibit a similar trend, with our 4x2 code outperforming both 
the DjABBA code and the BHV code at high SNR, and the DjABBA code in turn outperforming the 
BHV code. This can be attributed to the superior coding gain of our 4x2 code. The bad performance 
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of the BHV code at a high SNR is due mainly to the fact that it does not have full-diversity. Table 
Hm gives a comparison between the minimum determinants of the above three codes. The minimum 
determinants of our 4x2 code for 4-QAM and 16-QAM has been calculated using exhaustive search and 
the constellation energy hasn't been normalized to unity. However it has been ensured that the average 
energy per codeword has been maintained uniform for all the three codes. The DjABBA code that we 
have used for our simulations is the one that has been optimized for performance, and proposed in 
Chapter 9 of [10], . It can be seen that our code has a coding gain twice that of the DjABBA's. 

IX. Concluding Remarks 

In this paper, we have seen that it is possible to have full-rate codes with simplified ML-decoding 
complexity without having to sacrifice performance. We presented two codes, one each for the 2x2 
and the 4x2 MIMO system, both of which have lower ML-decoding complexity for general QAM 
constellations than the best known codes for such systems. Moreover, our 4x2 code outperforms the 
best DjABBA code while our 2x2 code matches the Golden code in performance. We also saw that 
the weight matrices play a decisive role in defining the ML-decoding complexity of an STBC and went 
on to show that some existing codes also offer simplified ML-decoding for square QAM constellations, 
something which was not known hitherto. Noting the similarity between the constructions of the 2x2 
code and the 4x2 code, it is natural to see if the design procedure can be extended to 2" transmit 
antennas, a > 2. However, there are two main issues to be concerned about: 

1) For our 2x2 code, we showed analytically that the minimum determinant for regular QAM 
constellations is 3.2. However, for our 4x2 code, we have checked that the minimum determinant 
for 4 and 16 QAM is 10.24 through exhaustive computer search. We couldn't do the same for 
higher constellation sizes, because such a search would run for weeks!. The rate of a square CIOD 
for 2°- transmit antennas is i^^, so that this STBC has 2a independent information symbols. If we 
were to extend our approach to 2°- transmit antennas, a > 2, the code would have 4a symbols and 
finding out the minimum determinant for 4 QAM itself would be time consuming. 

2) The ML-decoding complexity for our 2x2 code is of the order of M"^ and that for our 4x2 code 
is , for general constellations. So, the ML-decoding complexity for the STBC designed for 2°- 
transmit antennas, a > 2 would be of the order of M^"^^, while the rate would be While 
there is an increase in code rate, there is also a substantial increase in ML-decoding complexity, 
making the attractiveness of code design using this approach for higher number of transmit antennas 
questionable. 
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The following questions still remain unanswered. 

• For a 2 X 2 MIMO system, what is the minimum ML-decoding complexity achievable for a full- 
rate, full-diversity STBC ? Is it possible to have a full-rate, full-diversity code with an ML-decoding 
complexity of the order of Af ^ for all constellations. 

• Multi-group decodable codes [25] offer simplified ML-decoding complexity. For a given transmit 
antenna, what is the maximum rate that a multi-group decodable code can have ? For the 4x2 
MIMO case, is it possible to have a full-rate, full-diversity, two-group decodable STBC, so that the 
ML-decoding complexity is of the order of ? 
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Min det 


Code 


for M QAM 


Tilted QAM [20] 


0.8000 


Dayal-Varanasi code [8] 


3.2000 


The Golden code [7] 


3.2000 


Trace-orthonormal cyclotomic code [9] 


3.2000 


Paredes-Gershman code [11] 


2.2857 


Serdar-Sari code [12] 


2.0000 


The proposed code [22] 


3.2000 



TABLE I 

Comparison between the minimum determinants of some well known 2x2 STBCs 



Code 


ML Decoding complexity 


square QAM 


Non-rectangular QAM 


Tilted QAM 






Dayal-Varanasi code 


2A/VM 




The Golden code 






Trace-orthonormal cyclotomic code 






Paredes-Gershman code 




2M^ 


Serdar-Sari code 






The proposed code 


2M'^VM 


2Af3 



TABLE II 

Comparison between the ML-decoding complexity of some well known 2x2 STBCs for QAM 



Code 


Min det for 
4 and 16 QAM 


ML Decoding complexity 


Square QAM 


Non-rectangular QAM 


DjABBA code [10] 


0.64 


4M« 


2^ 


BHV code [17] 





4MVM 


2M^ 


The proposed code [24] 


10.24 




4M5 



TABLE III 

Comparison between the minimum determinant and the ML-decoding complexity of 4 x 2 STBCs for QAM 

constellations 
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Fig. 1. CER performance of 2 x 2 codes for 4-QAM 
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Fig. 2. CER performance of 2 x 2 codes for 16-QAM 
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Fig. 3. CER performance of 4 x 2 codes for 4-QAM 
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Fig. 4. CER performance for 4 x 2 codes for 16-QAM 
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32 QAM constellation 
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Fig. 5. 32 QAM - An example of a non-rectangular QAM constellations 
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